Short-range order of Br and three-dimensional magnetism in (CuBr)LaNb 2 07 
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We present a comprehensive study of the crystal structure, magnetic structure, and micro- 
scopic magnetic model of (CuBr)LaNb2C>7, the Br analog of the spin-gap quantum magnet 
(CuCl)LaNb2C>7. Despite similar crystal structures and spin lattices, the magnetic behavior and 
even peculiarities of the atomic arrangement in the CI and Br compounds are very different. The 
high-resolution x-ray and neutron data reveal a split position of Br atoms in (CuBr)LaNb207. This 
splitting originates from two possible configurations developed by [CuBr] zigzag ribbons. While the 
Br atoms are locally ordered in the ab plane, their arrangement along the c direction remains partially 
disordered. The predominant and energetically more favorable configuration features an additional 
doubling of the c lattice parameter that was not observed in (CuCl)LaNb207. (CuBr)LaNb207 
undergoes long-range antiferromagnetic ordering at Tjv = 32 K, which is nearly 70 % of the lead- 
ing exchange coupling J4 ~ 48 K. The Br compound does not show any experimental signatures 
of low-dimensional magnetism, because the underlying spin lattice is three-dimensional. The cou- 
pling along the c direction is comparable to the couplings in the ab plane, even though the shortest 
Cu-Cu distance along c (11.69 A) is three times larger than nearest-neighbor distances in the ab 
plane (3.55 A). The stripe antiferromagnetic long-range order featuring columns of parallel spins in 
the ab plane and antiparallel spins along c is verified experimentally and confirmed by the micro- 
scopic analysis. 



PACS numbers: 61.66.Fn, 75.30.Cr, 75.30.Et, 61.72.Ff 

I. INTRODUCTION 

Transition-metal halides are one of the best play- 
grounds for studying diverse magnetism of low- 
dimensional spin systems. Considering four stable halo- 
gen elements (F, CI, Br, and I), chlorine and bromine 
are most appealing because of the similar chemistry, 
yet different ionic radii that strongly influence superex- 
change pathways and relevant magnetic interactions. 
The change in the coupling regime between a chloride and 
an isostructural bromide is a fairly common scenario?"^ 
The differences between chlorides and bromides can be 
generally understood in terms of interatomic distances 
and angles, depending on the size of the ligand. While 
these geometrical parameters determine individual ex- 
change couplings, the influence on the spin lattice may 
be more involved, especially in mixed CI/Br systems 
with a variable crystallographic symmetry^ and effects 
of bond randomness P In the following, we will present 
an even more striking example, where Br atoms not 
only change the magnetic ground state but also trig- 
ger the incomplete structural order in the compound. 
We will consider (CuBr)LaNb207, a long-range ordered 
antiferromagnetJS which is remarkably different from its 
CI analog (CuCl)LaNb207 showing a gapped singlet 
ground state with zero ordered magnetic moment on the 
Cu sitesP 

Both (CuX)LaNb 2 7 compounds (X = CI, Br) feature 
flat [Cu +2 X] magnetic layers separated by non-magnetic 




FIG. 1. (Color online) Left panel: low-temperature crys- 
tal structure of (CuCl)LaNb207 featuring the [CuCl] zigzag 
ribbons and tilted Nb06 octahedral^ Right panel: spin lat- 
tice in the ab plane. Although the spatial arrangement of Cu 
atoms resembles the square lattice, numerous inequivalent in- 
teractions are present!^ Note the spin dimers formed on the 
fourth-neighbor Cu atoms, as shown by a shaded oval. For 
notation of individual couplings, see Fig. |14| 

[LaNbJ 5 07] slabs with two layers of corner-sharing NbOg 
octahedra (Fig. [I]). The original structural model de- 
scribed by Kodenkandath et a/P is deceptively simple, 
with a regular square-lattice arrangement of Cu and halo- 
gen atoms, and the ensuing tetragonal symmetry. Recent 
studietP^ questioned the completeness of this structural 
model, and identified an orthorhombic superstructure re- 
lated to the doubling of both the a and b lattice param- 
eters. In (CuCl)LaNb207, main features of this super- 
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structure are: i) cooperative tilts of the NbOg octahe- 
dra; and ii) displacements of Cu and CI atoms in the ab 
plane. The Jahn- Teller distortion inherent to Cu +2 re- 
sults in shorter and longer Cu-Cl bonds, so that Cu and 
CI atoms fo rm zigzag ribbons running along the b direc- 
tion (Fig. EESllillSI The orthorhombic structural model 
was derived from x-raj^ and neutron^ powder data, and 
later confirmed in a single-crystal x-ray experiment!^ 
This model satisfactorily accounts for the spin-gap mag- 
netic behavior of (CuCl)LaNb20y (note spin dimers in 
the right panel of Fig. [l]), altho ugh details of the spin 
lattice remain controversial liSilH 

In contrast to the CI compound, its Br analogue 
(CuBr)LaNb2 07 still lacks a detailed structural and 
microscopic investigation. Experimental data evidence 
magnetic couplings of J ~ 30 — 40 K, both ferromag- 
netic (FM) and antiferromagnetic (AFM)™ and a stripe 
(columnar) AFM long-range order that sets in below 
Tjv — 32 K. 6 While these observations would conform 
to the square-lattice arrangement of Cu +2 atoms with 
the dominating next-nearest-neighbor AFM coupling,™ 
an attempt to quantify the square-lattice model essen- 
tially failed and predicted a strongly frustrated regime 
of (CuBr)LaNb207, in agreement with naive computa- 
tional results™ and in apparent contradiction with the 
high Neel temperature of the compound (Tjy/J — 1)P 

Yoshida et aZ.™ used nuclear magnetic resonance 
(NMR) to demonstrate sizable deviations from the 
tetragonal symmetry at the Cu and Br sites, thus in- 
validating the available structural model as well as the 
square-lattice magnetic model. They also proposed a 
low-symmetry struct ure c ommon for (CuCl)LaNb 2 7 
and (CuBr)LaNb207j 17 * 18 l but recent diffraction studies 
of the CI compound refuted this conjectureP™ R en anc ] 
Chen™ p U t forward an alternative orthorhombic struc- 
ture later confirmed for (CuCl)LaNb207 and - as we will 
show below - roughly matching the experimental struc- 
ture of (CuBr)LaNb207, which however features a dis- 
ordered Br position (compare Figs. [I] and [5). No micro- 
scopic description for the magnetism of (CuBr)LaNb2 07 
has been given so far. 

In this paper, we first establish the crystal structure 
of (CuBr)LaNb 2 07 and further consider the microscopic 
magnetic model, thus providing a comprehensive descrip- 
tion of this compound in terms of the precise atomic 
positions and ensuing exchange couplings. We show 
that (CuBr)LaNb 2 07 is weakly frustrated and, more- 
over, reveals only weak quantum fluctuations because of 
the enhanced connectivity of the spin lattice. Our re- 
sults are naturally separated into the crystal-structure 
(Sec. III| and microscopic-model parts (Sec. IV) followed 
by a discussion and summary in Sec. |Vj where we com- 
pare (CuBr)LaNb207 to the isostructural compounds 
and consider broader implications for other Cu +2 bro- 
mides. 



II. METHODS 

Powder samples of (CuBr)LaNb 2 07 were prepared by 
a two-step procedure following Ref. [H First, we synthe- 
sized RbLaNb 2 07 by mixing stoichiometric amounts of 
La203 and M^Os with a 25% excess of Rb 2 C03, fol- 
lowed by an annealing at 850 °C (1123 K, 8 hours) and 
1050 °C (1323 K, 30 hours) with one intermediate re- 
grinding. The La 2 3 and Nb 2 5 powders were dried at 
800 °C prior to the experiment. The RbLaNb207 pow- 
der was washed with water to remove excess rubidium 
oxide. 

On the second step, RbLaNb 2 07 was mixed with a 
twofold amount of anhydrous CuBr 2 , pressed into a pel- 
let, sealed into an evacuated quartz tube, and heated at 
350 °C (623 K) for 48 hours. The resulting sample was 
again washed with water to eliminate the excess of CuBr 2 
as well as RbBr formed during the reaction. The qual- 
ity of starting materials, the RbLaNb207 precursor, and 
the final (CuBr)LaNb207 product was carefully checked 
with powder x-ray diffraction (XRD) measured using Hu- 
ber G670 Guinier camera (CuK Q i radiation, ImagePlate 
detector, 28 = 3 — 100 deg. angular range). The success- 
ful refinement of the high-resolution x-ray and neutron 
data (Sec. Ill I further confirms the high purity of our 
powder samples. 

High-resolution XRD data for structure refinement 
were collected at the ID31 beamline of European Syn- 
chrotron Radiation Facility (ESRF, Grenoble, France) 
using a constant wavelength of A ~ 0.4 A and eight 
scintillation detectors, each preceded by a Si (111) an- 
alyzer crystal, in the angular range 29 = 1 — 40 deg. The 
powder sample was contained in a thin-walled borosili- 
cate glass capillary that was spun during the experiment. 
The temperature of the sample was controlled by a He- 
flow cryostat (temperature range 10 — 200 K), a liquid- 
nitrogen cryostream (200 — 350 K), and hot-air blower 
(350 - 750 K). 

Neutron diffraction data were collected at the E9 high- 
resolution diffractometer (A ~ 1.797 A, Q = 1.2 — 
6.4 A" 1 , T = 20 K and 300 K) installed at the Hahn- 
Meitner Institute (Helmholtz Center Berlin for Materials 
and Energy, Germany). Unfortunately, these data were 
somewhat limited in the q range and insufficient to detect 
weak magnetic reflections. Therefore, we performed fur- 
ther experiments at the high-resolution instrument D2B 
(A ~ 1.595 A, Q = 0.9-7.6 A" 1 , T = 10 K) and the high- 
flux instrument D20 (A ~ 2.417 A, Q = 0.35 - 4.9 A" 1 , 
T = 1.5 K and 40 K), both installed at Institute Laue- 
Langevin (ILL, Grenoble, France). A 6 g powder sam- 
ple used for all neutron measurements was loaded into 
a cylindrical vanadium container and cooled down with 
the standard Orange He- flow cryostat. Crystal and mag- 
netic structures of (CuBr)LaNb2 07 were refined with 
JANA200(j™ anc j FullProi™ programs, respectively. 

The specimen for the electron diffraction (ED) study 
was prepared by crushing the sample under ethanol and 
depositing a drop of suspension on a holey carbon grid. 
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FIG. 2. (Color online) Left panel: high-resolution XRD patterns of (CuBr)LaNb207 measured at 20 K, 570 K, and 725 K for 
the a-, f3-, and 7-polymorphs, respectively. Only the superstructure reflections are indexed. The dashed line is the Rietveld 
refinement of the 20 K pattern. Note the diffuse scattering that is missing in the refinement. The diffuse features can be 
ascribed to a unit cell with the doubled c parameter. This diffuse scattering is denoted by arrows and persists at 570 K in the 
/3-polymorph, whereas the sharp superstructure reflections with even h + k disappear. Right panel: a part of the 20 K pattern 
showing the splitting of the 400 reflection. 



ED patterns were taken at room temperature with a 
Tecnai G2 transmission electron microscope operated at 
200 kV. Owing to the instability of the compound under 
intense electron beam, high-resolution imaging was not 
possible. 

Thermogravimetric analysis (TGA) was performed 
with the STA409 Netzsch thermal balance in the 300 - 
1050 K temperature range. For differential scanning 
calorimetry (DSC), we used a Perkin Elmer 8500 instru- 
ment in the temperature range 100 — 670 K. In both 
measurements, a heating rate of 10 K/min, corundum 
crucibles, and Ar atmosphere were used. 

The magnetic susceptibility of (CuBr)LaNb20y was 
measured with Quantum Design MPMS SQUID magne- 
tometer in the temperature range 2 — 380 K in applied 
fields between 0.1 T and 5 T. 

To evaluate lattice energies, individual exchange 
couplings, and the microscopic magnetic model, 
we performed full-potential scalar-relativistic density- 
functional-theory (DFT) band-structure calculations us- 
ing the FPL09 . 01-35 codeP^We used both the local den- 
sity approximation (LDA j^2l and the generalized gradient 
approximation (GGA^P^for the exchange-correlation po- 
tential. The symmetry-irreducible part of the first Bril- 
louin zone was sampled by a fine mesh of up to 570 points 
for the 48-atom crystallographic unit cell, 48 points for 
96-atom (doubled) supercells, and 8 points for 144-atom 
(tripled) supercells. Convergence with respect to the k 
mesh was carefully checked. Correlation effects in the 
Cu 3d shell were treated either on the model level (Hub- 
bard model constructed on top of LDA band structure) 
or in the framework of the mean-field DFT+C/ approach. 
Structure relaxations for large supercells were also per- 
formed in the VASP code 24 with the basis set of projected 
augmented waves^and the energy cutoff of 400 eV. Fur- 
ther details of the computational procedures are given in 



the respective sections. 

The microscopic magnetic model was further refined 
against the experimental data using quantum Monte 
Carlo (QMC) simulations of the magnetic susceptibil- 
ity, the magnetization isotherms, the Neel temperature, 
and of the ordered magnetic moment. These simulations 
were performed with the loorpH and dirloop_ss al- 
gorithms of the ALPS simulation package, 28 as further 
described in Sec. IIV CI 



III. CRYSTAL STRUCTURE 
A. Low-temperature structure 

Low-temperature x-ray and neutron diffraction pat- 
terns of (CuBr)LaNb 2 7 (Figs. [2] and [j§, as well as 
electron diffraction patterns measured at room temper- 
ature (Fig. [4jp2l showed weak reflections violating the 
tetragonal unit cell proposed in earlier studiesP- Addi- 
tionally, the high-resolution x-ray data revealed a weak 
orthorhombic distortion, evidenced by the splitting of the 
400 reflection shown in the right panel of Fig. [2j Most of 
the superstructure reflections could be indexed in a C- 
centered pseudotetragonal 2a su b x 2a su b x c su b unit cell, 
where a su b — 3.9 A and c su b — 11.7 A are parameters of 
the tetragonal subcell reported in Ref. [H The data also 
showed several diffuse features (Figs. [2] and [3]) along with 
faint electron-diffraction spots violating the C-centering. 
These features originate from the short-range order of Br 
atoms and will be discussed separately in Sec. |IIIB| 

Since no reflection conditions other than the general re- 
quirement h + k = 2n for the C-centering were observed, 
the Cmmm space group or its subgroups could be used 
for the structure refinement. In the Cmmm space group, 
the origin matches the inversion center that can be placed 
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FIG. 3. (Color online) Neutron diffraction patterns of 
(CuBr)LaNb207 measured at 1.7 K and 40 K, below and 
above the magnetic transition at Tjv — 32 K, respectively. 
Only the superstructure reflections are indexed. Note two 
weak nuclear reflections that can be indexed with the dou- 
bled c parameter, only. Magnetic reflections are labeled with 
arrows. The patterns are offset for clarity. 



TABLE I. Lattice parameters (in A) and space groups for 
structure refinements of (CuBr)LaNb207 at different temper- 
atures T (in K), as given in Tables [Til and [!vl The error bars 
are based on the Rietveld refinement. 



T 


x/a 


y/b 


z/c 


Space group 


20 


7.7856(2) 


7.7983(2) 


11.6938(2) 


Cmmm 


300 


7.7926(2) 


7.8009(2) 


11.7017(2) 


Cmmm 


720 


3.9077(1) 


3.9077(1) 


11.7192(1) 


PA/mmm 



either on Cu atoms, on Br atoms, or in between the two 
neighboring Cu atoms. Only the latter setting allowed 
the successful refinement. In this case, the NbC>6 octa- 
hedra follow the a a b~c° tilting distortion,^ whereas the 
Cu and Br atoms are constrained to the 2 = plane by 
the respective mirror symmetry. Additionally, the m y 
and m x planes allow the displacements of the Cu and Br 
atoms along the a and b directions, respectively. While 
the a and b axes are interchangeable, we use the above 
setting with Cu on the m y plane and Br on the m x plane 
for the sake of consistency with the previously published 
data on (CuCl)LaNb 2 7 P 

The refinement of the high-resolution neutron data 
was somewhat unstable, because the weak orthorhombic 
splitting caused large fluctuations of the a and b lattice 
parameters. To improve the refinement, we additionally 
considered the high-resolution x-ray data, where the or- 
thorhombic splitting is better resolved (Fig. [2| . We have 
also used the neutron data from the high-flux D20 in- 
strument that is most sensitive to weak superstructure 
reflections. Altogether, the low-temperature structure 
of (CuBr)LaNb207 was simultaneously refined against 
three datasetsPH Diffuse scattering and weak magnetic 
reflections observed in the D20 experiment were excluded 



FIG. 4. Room-temperature electron diffraction patterns of 
(CuBr)LaNb207. Bright dots are the subcell reflections. Less 
bright dots with odd h and k are visible in the [001] and 
[1 10] patterns, and evidence the superstructure with the C- 
centered 2a su b x 2a su b x c su b unit cell. Faint dots with odd 
h + k violate the C-centered unit cell and originate from the 
short-range order of Br atoms (see text for details). 



from the refinement and analyzed separately (Sees. Ill B 
andlrVBl. 



The refinement of atomic positions and atomic dis- 
placement parameters (ADPs) produced an unusually 
large ADP for Br, Ui SO ~ 0.050 A 2 indicating a disorder of 
the Br atoms. A subsequent refinement of the anisotropic 
ADP showed a drastic elongation of the ellipsoid along 
the a direction. Therefore, a split position with the Br 
atom shifted from the m x plane was introduced (Fig. [5]). 
This improved the refinement and reduced the ADP of 
Br down to f.O x 10~ 2 A 2 . A displacement out of the 
m z plane allowed for a further reduction in the ADP and 
resulted in a non-zero z coordinate. Therefore, the Br 
atoms occupy the split 16-fold position in the averaged 
C-centered orthorhombic crystal structure. 

Final atomic positions are summarized in Table|Tl) 31 l 32 l 
Figure [5] depicts the low-temperature structure of 
(CuBr)LaNb2 07, which is further referred to as the a- 
polymorph. This structure bears apparent similarities 
to the low-temperature a-(CuCl)LaNb207 structure that 
also has the orthorhombic symmetry. Both CI and Br 
compounds reveal the a°b~c° tilting distortion with the 
out-of-phase rotation of the NbOg octahedra about the 
b axis. The tilting angle measured as the tilt of the 
Nb-03 bond with respect to the c axis is 4.8 deg in 
a-(CuBr)LaNb207 compared to 5.7 deg in the CI com- 
pound (Ref. llip . Following this rotation, terminal oxygen 
atoms displace along the a direction and promote similar 



TABLE II. Atomic positions and atomic displacement pa- 
rameters (U iso , in 10~ 2 A 2 ) for a-(CuBr)LaNb 2 7 at 20 K 
(first line) and 300 K (second line), as refined from the syn- 
chrotron XRD and neutron data. The space group is Cmmm. 
The Uiso of oxygen atoms were refined as a single parameter. 
The error bars are based on the Rietveld refinement. See text 



and Supplementary information for details. 


Atom 


Position 


x/a 


y/b 


z/c 


U iBO 


Cu 


4p 


0.7296(2) 


■i 

2 





0.36(3) 






0.7265(5) 


1 

2 





0.97(7) 


t-, □ 


16r 


0.4536(4) 


0.2385(3) 


0.0068(9) 


0.54(13) 






0.4568(8) 


0.2424(5) 


0.012(1) 


1.2(3) 


La 


4j 





0.2588(2) 


2 


0.02(2) 









0.2568(3) 


1 

2 


0.51(3) 


Nb 


So 


0.7467(2) 


2 


0.31005(4) 


0.00(2) 






0.7476(2) 


2 


0.30929(8) 


0.51(3) 


01 


4/ 





1 

2 


0.3265(4) 


0.20(1) 









1 
2 


0.3280(11) 


1.47(5) 


02 


8m 


i 

4 


3 
4 


0.3456(3) 


0.20(1) 






1 

4 


3 
4 


0.3428(9) 


1.47(5) 


03 


8o 


0.2276(2) 





0.15958(7) 


0.20(1) 






0.2352(7) 





0.1583(2) 


1.47(5) 


04 


4/i 


0.7831(3) 


1 

2 


i 

2 


0.20(1) 






0.7676(10) 


1 
2 


1 

2 


1.47(5) 


05 


4k 








0.3621(4) 


0.20(1) 












0.3634(10) 


1.47(5) 



Fixed occupancy g^ T 



displacements of Cu atoms in the ab plane (Fig. [5]) . 

An important difference between a-(CuCl)LaNb2 07 
and a-(CuBr)LaNb20y is the arrangement of halogen 
atoms. While the CI atoms are fully ordered and lie in the 
ab plane, the Br atoms are disordered in the ab plane and 
additionally show marginal out-of-plane displacements. 
This increases the symmetry from Pbam in the CI com- 
pound to Cmmm in the Br compound. Both structures 
reveal shortened Cu-X distances of about 2.39 A (X = 
CljPand 2.49-2.55 A (X = Br), compared to 2.76 A in 
the tetragonal structures reported earlierPThe elongated 
Cu-Br distances of 2.95 — 3.10 A are also present. The 
shortened and elongated bonds reflect the Jahn- Teller 
distortion driven by the d 9 electronic configuration of 
Cu+ 2 (Ref.[12]). 

Each Cu atom forms two short bonds to the 03 atoms 
of the NbC>6 octahedra (see Table [V]). Therefore, only 
two short bonds to Br may be formed to avoid over- 
bonding. In a-(CuCl)LaNb207, the irons-arrangement 
of the short Cu-Cl bonds (Cl-Cu-Cl angle close to 180°) 
is well defined, whereas the split Br position could allow 
for different local environments of Cu. However, elec- 
tronic effects strongly disfavor a random arrangement of 
Br and lead to the local ordering that largely resembles 
the cv-(CuCl)LaNb207 structure. Details of the short- 
range order in a-(CuBr)LaNb207 are discussed below 
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FIG. 5. (Color online) Low-temperature crystal structure 
of (CuBr)LaNb207 with the split position of Br atoms (the 
displacements along the c direction are too small to be visual- 
ized on this scale). Arrows show the displacements of oxygen 
atoms upon the a°b~c° tilting distortion (right panel) and 
the ensuing displacements of Cu atoms in the ab plane (left 
panel). The lines denote short Cu-Br bonds in the ab plane. 



in S ec. |IIIB| followed by microscopic considerations in 
Sec.HTDl 



B. Short-range order 

To explore the short-range structural order in a- 
(CuBr)LaNb207, we consider the NMR data reported 
by Yoshida et aZ.fEO who evaluated electric field gradients 
(EFGs) on Cu and Br sites. Their results evidence sin- 
gle crystallographic positions for both Cu and Br within 
the experimental resolution. The refined crystal struc- 
ture allows for several local configurations (Fig. |6J): i) 
irons-arrangement of short Cu-Br bonds forming zigzag 
ribbons along the b direction; ii) cis-arrangement of short 
bonds forming dimers (similar to the earlier proposal 
in Ref. 1 1T[) : iii) cis-arrangement of short bonds forming 
chains along b. The EFGs calculated for these config- 
urations along with the experimental results are listed 
in Table [TTTJ In DFT calculations, we use the LSDA+C/ 
method with the on-site Coulomb repulsion Ud = 5 eV, 
the on-site Hund's exchange Jd = 1 eV, and the around- 
mean-field (AMF) double-counting correction. The vari- 
ation of these parameters in a reasonable range has little 
influence on the computed EFG values. Further details 
of the computational method are given in Sec. 



The EFG tensor is quantified by its leading component 
V zz and the asymmetry parameter r\ = (V yy — V xx )/V zz . 
Considering experimental accuracy and resolution, com- 
putational estimates show rather small differences be- 
tween three possible configurations. Although V zz of Cu 
changes sign depending on the cis- or irons-arrangement, 
this effect can not be captured experimentally, because 
NMR does not resolve the sign of V zz unambiguously. 
Nevertheless, experiments safely establish that only one 
local arrangement is present in the crystal structure, be- 
cause the respective V zz values on Cu sites differ for about 
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TABLE III. Calculated and experimental parameters of the 
EFG tensor for the Cu and Br sites in a-(CuBr)LaNb207. 
Vzz is the leading component (in 10 21 V/m 2 ) and r/ = (V yy — 
V xx )/V zz is the asymmetry. The models of local order are 
shown in Fig. [6] Model III features two inequivalent Cu sites 
with slightly different Cu-Br bond lengths. Experimental val- 
ues are taken from Ref. [T7| Note that the sign of V zz can not 
be determined experimentally. 



Model I (trans) 
Model II (cis) 
Model III (cis) 
Experiment 



Cu 

V zz r, 
-10.3 0.27 
12.4 0.34 
12.4/12.0 0.06/0.38 
±10.6 



Br 

v„ n 

-45.0 0.28 

-42.3 0.26 

-42.3 0.04 

±30.0 0.33 





Model I 



ribbon A 



ribbon B 





Model II 



Model III 



FIG. 6. (Color online) Possible local order of Br atoms in 
a-(CuBr)LaNb207. The top panel shows the only stable con- 
figuration (trans-arrangement of short Cu-Br bonds) and the 
principal scheme of disorder by random shifts of zigzag [CuBr] 
ribbons for one half of the b lattice parameter (compare to the 
refined structure shown in Fig. [5b. The b ottom panel shows 
the unstable configurations (see Sec. IIIBl. 



20 % and should be easily resolvable. According to the 
Cu V zz , the trons-arrangement of the Cu-Br bonds is 
more likely, as further confirmed by energetic consider- 
ations in Sec. |III D| All three configurations systemati- 
cally overestimate V zz for the Br site. Since a similar 
30 % overestimate has been reported for the V zz of CI 
in a-(CuCl)LaNb 2 7 (see Table V in Ref. the dis- 
crepancy is likely related to drawbacks of the computa- 
tional method and not to the inaccuracies of our struc- 
tural model. 



The trons-arrangement of short Cu-Br bonds results in 
the formation of planar CuC<2Br2 plaquettes and [CuBr] 
zigzag ribbons, similar to the [CuCl] ribbons in the struc- 
ture of a-(CuCl)LaNb 2 7 (see Fig. [TJ. The disorder in 
the Br position implies that such ribbons randomly dis- 
place for one half of the 6 lattice parameter to form an- 
other ribbon with a different arrangement of Br atoms 
(see the top panel of Fig. [6]). However, the ribbons 
within the ab plane should be ordered, because the shift 
of the ribbon for 6/2 creates Cu atoms with the cis- 
configuration of short bonds at the interface (see Fig.[l0|. 
Since the cis-arrangement of short bonds is energetically 
highly unfavorable (see Sec. 



NMR, the disorder within t 



III D I and is not observed by 



re ab plane can be ruled out. 



Thus, the only plausible scenario is the random arrange- 
ment of the [CuBr] ribbons in neighboring planes. 

Labeling one possible arrangement of the [CuBr] rib- 
bon as A and the ribbon displaced for 6/2 as B (Fig. [6j), 
one can construct structures with different orderings 
along the c direction. The simplest examples are AAAA 
(or the equivalent BBBB) and ABAB shown in Fig. [9] 
The former is equivalent to the a-(CuCl)LaNb20y struc- 
ture and has the 2a su b X 2a su b x c su b unit cell with the 
Pbam symmetry. The latter features the doubled trans- 
lation along c (2a su b x 2a su b x 2c su b cell). Ribbon A 
is transformed into ribbon B by a [^,§,0] translation 
supplied with the [0, 0, |] translation between the neigh- 
boring layers. The structure is, therefore, body-centered 
and retains all symmetry elements of the AAAA config- 
uration, thus resulting in the Ibam space group. While 
the AAAA structure should manifest itself by additional 
superstructure reflections violating the C-centering (i.e., 
hkl with odd h + k and integer I), the ABAB structure 
gives rise to the reflections with odd h+k and half-integer 
I. A closer inspection of the diffraction data shows that 
the reflections of the latter type are indeed present. 

The low-temperature neutron diffraction patterns 
(Fig. [3| feature two weak and somewhat diffuse reflec- 
tions that can be indexed with odd h + k and half- integer 
I. The missing intensity in the refinement of the low- 
temperature XRD data (Fig. [2| can be explained in a 
similar way. Finally, the faint spots at hkO (h + k = odd) 
in the [001] ED pattern (Fig. Eb are traces of the dif- 
fuse intensity lines intersecting the [001] reciprocal lat- 
tice plane. The x-ray data additional show very broad 
features at the positions of hkl reflections with odd h + k 
and integer I (e.g., the broad hump around 26 ~ 10.5 deg 
in Fig. [2]). This intensity may signal traces of the AAAA 
type (Pbam) ordering, which however is not seen in the 
neutron data. Since neutrons probe the bulk of the sam- 
ple, we conclude that in a-(CuBr)LaNb20 7 the ABAB 
type of the local order prevails, although sporadic inter- 
ruptions of this configuration give rise to the disordered 
position of Br in the averaged crystal structure with the 
Cmmm symmetry. Small regions of the AAAA type or- 
der that manifests itself in the x-ray and electron diffrac- 
tion data may also be present in the sample, although 
their concentration is low. Further studies on single crys- 
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tals should better resolve the diffuse intensity and provide 
additional insight into the short-range order of Br atoms 
in a-(CuBr)LaNb 2 7 . 

A refinement of the data in the Ibam space group ac- 
counts well for the intensities of all superstructure re- 
flections. However, it also shows the sizable broadening 
of the reflections with odd h + k, and arrives at a rela- 
tively large ADP of Br (U iso ~ 0.02 A 2 ). This confirms 
that the Br atoms are only locally ordered, and the aver- 
aged structure should be described by the Cmmra space 
group. The disorder in the positions of Br atoms along 
the c direction and the preference for the ABAB-type or- 
der are further confirmed by microscopic considerations 
in Sec. ImDl 



C. Temperature evolution 

In (CuCl)LaNb 2 07, the Cl-containing counterpart of 
(CuBr)LaNb207, the study of the temperature evolution 
provided a valuable insight into the formation of the su- 
perstructure, and clarified the interplay of displacements 
observed for different atomic speciesP Upon heating, the 
tilting distortion is eliminated above T\ ~ 500 K (a — > j3 
transition). Since the terminal oxygen atoms of the NbOe 
octahedra are no longer displaced, the ordered arrange- 
ment of Cu atoms is also lost (see Fig. [5]). However, 
the ordered positions of CI atoms are preserved up to 
the second — > 7 structural transition at T 2 ~ 640 K. 
Above T 2 , (CuCl)LaNb 2 07 displays the tetragonal crys- 
tal structure (a su b X a su b x c su b unit cell, 7-phase) with 
random displacements of both Cu and CI atoms within 
the [CuCl] planes. Note, however, that neither Cu nor CI 
occupy the high-symmetry positions on the four-fold axes 
and rather remain displaced, in order to retain the short 
Cu-Cl bonds (Jahn- Teller distortion) at least locallyP 

On the experimental side, the two-step formation of 
the a-(CuCl)LaNb 2 07 superstructure implies the ap- 
pearance of the hkl superstructure reflections in two 
stages. The superstructure reflections with odd h + k 
are observed in the /3-polymorph below T 2 (ordering of 
CI atoms), while the reflections with even h + k appear 
in the a-polymorph below T\ only (tilting distortion and 
ordering of Cu atoms). 

The temperature evolution of (CuBr)LaNb 2 07 follows 
the same line. The superstructure reflections with even 
h + k disappear above T\ ~ 500 K, whereas the diffuse 
scattering centered at hkl with odd h + k and non- integer 
I persists up to T 2 ~ 620 K (see Fig. [2). Although the 
Br atoms are locally ordered, their arrangement is not 
influenced by the onset of the tilting distortion in the 
[LaNb 2 C>7] perovskite slabs. 

Upon heating, the weak orthorhombic distortion 
present in (CuBr)LaNb 2 07 is gradually reduced and 
eventually disappears around T\ (Fig. [7]). Unfortunately, 
this effect can not be tracked precisely because of the 
sizable reflection broadening that prevented us from the 
independent refinement of the a and b parameters above 



TABLE IV. Atomic positions and atomic displacement pa- 
rameters (Uiso, in 10 -2 A 2 ) for 7-(CuBr)LaNb207 refined us- 
ing the synchrotron XRD data at 720 K. The space group is 
Pi/mmm. The Ui so of oxygen atoms were refined as a single 
parameter. The error bars are based on the Rietveld refine- 
ment. See text and Supplementary information for details. 
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lb 
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FIG. 7. (Color online) Temperature evolution of subcell 
parameters and subcell volume for (CuBr)LaNb20 7 . Lines 
are guide for the eye. 



450 K. Above T 2 , the c parameter becomes temperature- 
independent, whereas the temperature evolution of the 
cell volume shows a slight change in the slope. DSC does 
not reveal any anomalies at T± and T 2 , because the struc- 
tural changes are weak and associated with only marginal 
changes in the entropy. TGA confirms the stability of 
(CuBr)LaNb 2 07 within the studied temperature range 
and shows the onset of decomposition (weight loss) above 
750 K, onlyP 

Above T 2 , the crystal structure of (CuBr)LaNb 2 07 can 
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FIG. 8. (Color online) Temperature-dependent out-of-plane 
displacement of the Br atoms (AzBr), as obtained from the 
refinement of the high-resolution XRD data and the joint re- 
finement of x-ray and neutron data (XRD+NPD). The inset 
shows the [CuBr] ribbons with the split position of Br above 
and below the mirror plane at z — 0. 



be refined in a tetragonal a su b X a su b X c su b unit cell 
(7-polymorph) The regular [LaNb^Oy] slabs have the 
four-fold symmetry, while the Cu atoms are disordered 
over four equivalent positions in the z = plane, and the 
Br atoms are disordered over eight equi vale nt positions 
lying above and below this plane (Table [iv]). 31 The dis- 
placements of Cu and Br atoms ensure the formation of 
short Cu-Br bonds that preserve local Jahn- Teller dis- 
tortions for Cu atoms. The out-of-plane displacement 
of Br [AzBr = zbt/c = 0.0189(5)] is notably increased 
compared to Aze r — 0.007(1) at low temperatures. This 
feature is clearly different from 7-(CuCl)LaNb 2 07, where 
CI atoms remain in the mirror plane and have a larger 
ADP (U iso = 0.047(2) A 2 ) with a stronger tendency to 
the in-plane displacements 

The out-of-plane Br displacements systematically in- 
crease with temperature (Fig. [8]), as confirmed by Ri- 
etveld refinements of the XRD data and the joint re- 
finement of the x-ray and neutron data at room temper- 
ature (Tabic |n| . The three- fold increase in Azer upon 
heating shows that the displacements are largely driven 
by dynamic effects, such as soft phonon modes. How- 
ever, at low temperatures the sizable displacements are 
still present and manifest one of the differences between 
(CuBr)LaNb207 and its Cl-containing counterpart. 

The smooth temperature evolution of lattice parame- 



ters and cell volume (Fig. |III C ) suggests that the phase 
transitions at T% and T 2 are of second order, similar to 
(CuCl)LaNb 2 7 (Ref.E§. 



D. Energetics 

We will now consider the crystal structure of 
(CuBr)LaNb2 07 from a microscopic viewpoint, and com- 
pare total energies for different arrangements of Br 
atoms. First, we investigate two ordered models (Fig. [9]) 




Pbam (AAAA) 



Ibam (ABAB) 



FIG. 9. (Color online) Structural models representing the 
short-range order of Br atoms in (CuBr)LaNb207. The pre- 
vailing Ibam model (right panel) results in the doubling of the 
c lattice parameter and manifests itself by diffuse scattering 
in the x-ray, neutron, and electron diffraction data. 



that are consistent with the experimental diffraction data 
reported in Sections |III A| and |IIIB| One model repre- 
sents the Pbam structure of cv-(CuCl)LaNb 2 07 type with 
the c su b lattice parameter and the equivalent arrange- 
ment of [CuBr] ribbons along the c direction (AAAA). 
The other model features the Ibam symmetry and the 
doubled lattice parameter along c because of the alter- 
nating arrangement of the [CuBr] ribbons (ABAB). Both 
models were optimized^ until residual forces dropped 
below 0.01 eV/A. The ensuing interatomic distances and 
angles are listed in Table [V] along with the experimen- 
tal results for the low-temperature crystal structure from 
Table IE 

The optimized crystal structures show remarkable 
agreement with the experimental data, thus suggesting 
an excellent potential of DFT+C7 for structure predic- 
tion. The different arrangement of the [CuBr] ribbons in 
the Pbam and Ibam models has no appreciable influence 
on individual interatomic distances and angles. Indeed, 
both structures have similar energy, with a slight prefer- 
ence for the Ibam model (~ 2 meV/f.u.) in agreement 
with the experimental data (Sec. |IIIB"| ). The marginal en- 
ergy difference between the Pbam and Ibam structures is 
a natural explanation for the incomplete structural order 
in (CuBr)LaNb 2 7 . 

For the sake of completeness, we also estimated the en- 
ergies for other possible configurations of Br atoms. The 
structures with the czs-arrangement of short bonds, as 
depicted in Fig. [6j appeared to be highly unstable and 
converged to the structure with the irans-arrangement. 
However, the cis-arrangement of short bonds can be sta- 
bilized as a defect in the regular structure with Cu atoms 
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TABLE V. Experimental (at T = 20 K) and relaxed 
interatomic distances (in A) and angles (in deg) in the 
(CuBr)LaNb207 structure. The structures were relaxed us- 
ing GGA+[/ with Ud = 5 eV. Note that the experimental 
structure features a split position of Br, hence the number of 
Cu-Br distances is doubled. 





Experiment 


Relaxation 


Relaxation 


Space group 


Cmmm 


Pbam 


Ibam 


Cu-03 


1.866(1) 


1.848 


1.848 


Cu-Br 


2 x 2.490(3) 


2.517 


2.516 


Cu-Br 


2 x 2.550(3) 


2.545 


2.543 


Nb-Ol 


1.981(2) 


1.968 


1.968 


Nb-02 


1.994(1) 


2.000 


2.005 


Nb-02 


1.994(1) 


2.005 


2.001 


Nb-03 


1.766(1) 


1.788 


1.788 


Nb-04 


2.239(1) 


2.271 


2.271 


Nb-05 


2.015(2) 


2.046 


2.046 


Cu-Br-Cu 


101.8(2) 


101.6 


101.7 



featuring a more stable trans-arrangement. To intro- 
duce such a defect, we used the unit cell doubled or 
tripled along the a direction, and created inequivalent 
[CuBr] ribbons. Fig. 10 depicts the fully relaxed struc- 
ture with the cis-arrangement of Cu-Br bonds on the 
interface. Although such structures are metastable and 
form local minima on the potential energy surface, they 
are highly unfavorable with respect to the Pbam and 
Ibam structures shown in Fig. [9j The relaxations per- 
formed for the doubled and tripled unit cells show ex- 
cess energies of 0.158 eV/f.u. and 0.125 eV/f.u., respec- 
tively, thus yielding the energy cost of 8 x 0.158/4 eV 
and 12 x 0.125/4 eV or 0.30 - 0.35 eV for each Cu atom 
with the czs-arrangement of short bonds. This energy 
cost is significantly higher than the thermal energy of 
0.055 eV available at the preparation temperature of 
600 — 650 K. Therefore, the formation of Cu atoms with 
the cis-configuration of short bonds is extremely unlikely. 

Our structure relaxations show that different arrange- 
ments of Br atoms have little influence on the tilting 
distortion developed by the [LaNb 2 07] perovskite slabs. 
The formation of defects in the [CuBr] layers (Fig 



10) 



does not disrupt the uniform a &~c u tilting pattern. This 
illustrates the remarkable independence of the Br dis- 
placements and the tilting distortion that set in upon 
two different phase transitions and remain essentially de- 
coupled. 

Finally, we examined the possible displacements of Br 
atoms out of the z = plane. Structures with the 
out-of-plane displacements lack the m z mirror plane and 
have either P2 (P112) or PI symmetries. Both space 
groups feature two inequivalent Br positions, Br a and 
Brh, and allow for two different out-of-plane displace- 
ments each. The structures with Br a and Brj, displaced in 
the same/opposite directions with respect to the z — 
plane are denoted ft an d t4-j respectively. Altogether, 




FIG. 10. (Color online) Sketch of the possible disorder in 
the ab plane: the A- and B-type [CuBr] ribbons coexist within 
the same layer of the (CuBr)LaNb207 structure. Light and 
dark lines show the short Cu-Br bonds in the trans- and cis- 
configurations, respectively. Note that the change in the con- 
figuration of the ribbons does not disrupt the tilting distortion 
in the [LaNb2 0r] slabs (bottom panel). 



four different starting configurations are possible (see the 
bottom panel of Fig. [IT] ), but all of them relax to the 
same Pbam (or Ibam, depending on the chosen order- 
ing type) structure with zb? = and the mirror-plane 
symmetry restored. The lowest energy of the Pbam and 
Ibam structures can also be seen from the upper panel of 
Fig. [IT] where total energies for fully relaxed structures 
with fixed Br displacements are presented.^ 

The out-of-plane Br displacements are not equivalent 
in terms of the energy (Fig. 11). The displacements re- 
quire lower energies in the PI ft an d P2 fi structures, 
where the two Br atoms of the same Cu02Br2 plaque- 
tte are shifted in the opposite directions. This way, the 
whole plaquette is slightly tilted, while its overall shape 
is preserved. The displacements in the PI fl an d P2 
ft structures require much higher energies, because the 
shape of the plaquette is changed. Finally, the preference 
of the PI ft structure compared to the P2 structure 
can be explained by the increase in the Br-Br distances 
between the negatively charged Br ions of neighboring 
[CuBrl ribbons. 



The out-of-plane displacements require relatively low 
energies even for the unfavorable configurations PI 
and P2 ft- For example, the thermal energy of about 
40 meV at 500 K is sufficient to induce random displace- 
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FIG. 11. (Color online) Top panel: relative energies for dif- 
ferent out-of-plane displacements of Br atoms, with all other 
atomic positions fully optimized. Bottom panel: sketch of 
the Br displacements, see text for notation details. Different 
colors show the Br a and Br;, sites, as explained in the text. 



merits with Aze r — ±0.02. This estimate compares well 
to the experimental displacement of Aze r — 0.02 ob- 
served at high temperatures (Fig. [8|. Our microscopic 
study confirms the dynamic nature of the out-of-plane 
Br displacements. At low temperatures, the displace- 
ments should be eliminated to restore the mirror-plane 
symmetry of the equilibrium crystal structure. The resid- 
ual out-of-plane displacements observed experimentally 
at low temperatures are a non-equilibrium feature, pre- 
sumably related to defects and/or incomplete ordering of 
Br atoms. 



IV. MICROSCOPIC MAGNETIC MODEL 
A. Electronic structure 

To evaluate the microscopic magnetic model, we calcu- 
late the band structure and evaluate individual exchange 
couplings. Our results are based on the fully relaxed 
crystal structures with Pbam or Ibam symmetry (Fig. [9]) . 
Since interatomic distances and angles in these structures 
are nearly indistinguishable (see Table [v} , the ensuing 
exchange couplings are also the same. Reference calcula- 
tions for ordered structures based on the experimental re- 
finement (neglecting the out-of-plane displacements of Br 
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12. (Color online) LDA density of states for the low- 



temperature structure of (CuBr)LaNh>207. 
is at zero energy. 



The Fermi level 



atoms) produced similar results for individual exchange 
couplings. 

The LDA band structures of (CuBr)LaNb2 07 and 
(CuCl)LaNb207 are rather similar (compare Fig. 12 to 
Fig. 2 of Ref. IT4|) . Oxygen 2p bands span the energy 
range between —7 and —1 eV, with Cu 3d bands above 
—4 eV and Br Ap bands above —1.5 eV. The states at 
the Fermi level are represented by four narrow bands ac- 
cording to four Cu atoms in the unit cell. These bands 
are formed by Cu d x i_ y i orbitals with an appreciable 
admixture of the Br Ap and O 2p orbitals (here, x and 
y are directed along the Cu-0 and short Cu-Br bonds, 
respectively) . The sizable energy separation between Cu 
d x 2_ y 2 and the rest of Cu 3d states is driven by the large 
crystal-field splitting that illustrates the Jahn- Teller dis- 
tortion underlying the formation of short Cu-Br bonds. 
According to Ref. Q21 the simple tetragonal crystal struc- 
ture from Ref. |5| entails Cu atoms with four equivalent 
Cu-Br bonds. The respective orbital state is nearly de- 
generate, because the two e g orbitals have similar ener- 
gies in the effectively octahedral Cu02Br4 coordination. 
The displacements of Br atoms off the four-fold rota- 
tion axis and the ensuing shortening of the two Cu-Br 
bonds stabilize the non-degenerate orbital configuration. 
In (CuBr)LaNb 2 07, the contribution of halogen p states 
at the Fermi level is slightly larger than in the CI com- 
pound, 36 % and 33 %, respectively. This trend follows 
the higher energy of Br Ap compared to the CI 3p orbitals. 

The dispersions of the Cu d x 2_ y 2 bands identify pos- 
sible electron transfers and magnitudes of AFM interac- 
tions in the system (Fig. 13 1. The electron transfers (U) 



are quantified in a tight-binding model based on Wannier 
functions with the d x 2_ y 2 orbital character We further 
supply this tight-binging model with a Hubbard term to 
account for an effective on-site Coulomb repulsion U e g, 
and reduce the problem to the strongly correlated limit 
tj -C U c ff. This way, one arrives at the Heisenberg model 
with AFM exchanges jf" FM = 4tf/U e g for low-lying ex- 
citations. We use U c h = 4 eV, according to previous 
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FIG. 13. (Color online) LDA band structure in the vicin- 
ity of the Fermi level (thin light lines) and the fit with 
the tight-binding model (thick dark lines). The Fermi level 
is at zero energy. The notation of k points is as fol- 
lows: r(o,o,o), 0,0), M(|,i,0), F(0,|,0), Z(0,0,|), 

M2>0j 2)' ^(2 J 21 2")' ^-(0' 2> 2")' 



TABLE VI. Magnetic couplings in (CuBr)LaNb 2 7 : Cu- 
Cu distances, hopping parameters i, of the tight-binding 
model together with the ensuing AFM exchange couplings 
jAfm _ 4^2 /f/ cff 5 anc i f u n exchange couplings Ji evaluated 
using DFT+C calculations with the AMF and FLL double- 
counting correction schemes. See text for details. 
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Ji (AMF) 


Ji (FLL) 
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(eV) 
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Jl 
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-75 


-47 


Ji 


3.575 


0.045 


24 


31 
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4.211 


0.034 


13 
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■h 


5.510 


-0.047 


26 
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5.510 


0.011 


1 


-4 


-14 


Js 


7.786 


-0.004 









Jk 


7.798 


0.008 


1 






Ja 


8.579 


-0.097 


110 


144 


54 


J'a 


8.862 


-0.036 


15 


37 


17 


■h 


11.020 


0.011 


1 






Jk 


11.020 


0.035 


14 


13 a 


1° 


J± 


11.694 


-0.028 


9 


17 
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Here, only the sum J5 + JL was evaluated. 



studies of Cu +2 halides and oxyhalides! 14 * 38 ! 

The hoppings of the tight-binding model and the re- 
sulting AFM exchange couplings are listed in Table |VI| 
Following Ref. [21 we denote the couplings in the ab 
plane according to the Cu-Cu distances (Ji, J[, and J" 
between nearest neighbors, J 2 and J' 2 between second 
neighbors, etc.) and use Jj_ for the interlayer coupling 
along c (Fig. 14 1. At first glance, the microscopic sce- 
nario is somewhat unusual, because the leading coupling 
J 4 runs between fourth neighbors, while the nearest- 
neighbor and next-nearest-neighbor couplings are rela- 
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FIG. 14. (Color online) Magnetic couplings in the ab 
plane (left panel) and the respective spin lattice (right panel). 
Light and dark spheres denote up and down spins in the or- 
dered antiferromagnetic state, as determined experimentally 
(Sec. IV B I. The couplings J2, J5, and J5 are not shown. The 
interlayer coupling J± runs along the c direction. 



tively weak. This peculiar behavior is rationalized in 
terms of the magnetic d x 2_ y 2 orbital lying in the plane 
of the Cu02Br2 plaquette. The large spatial extent of 
Br Ap orbitals facilitates the long-range hopping between 
the fourth neighbors, but this effect is highly selective. 
Only the fourth-neighbor hoppings along the short Cu- 
Br bonds are enhanced (e.g., along [210] and not along 
[2l0], or the other way around, depending on the orien- 
tation of the plaquette). Moreover, the t± path is twice 
more efficient than its t\ counterpart. This can be ex- 
plained by the different Cu-Br-Br angles of 156.4 deg 
and 144.7 deg, respectively. The larger angle renders 
a more straight and, therefore, more efficient superex- 
change pathway (see the left panel of Fig. 14). Note that 
the interlayer hopping t± is comparable to several short- 
range intralayer interactions, such as t[, t", t 2 , and t\. 
Similar to the Cu-Br-Br-Cu pathway for J4, the long 
Cu-O-Nb-O-Nb-O-Cu pathway is also rather efficient 
despite the very long Cu-Cu distance of 11.69 A, be- 
cause oxygen 2p orbitals entering the Wannier functions 
are directed along this pathway, while the low-lying Nb 
Ad states additionally facilitate the superexchange. 

Altogether, the LDA-based microscopic model of 
(CuBr)LaNb2 07 strongly resembles the respective 
model for the isostructural CI compound, where 
the leading coupling between fourth neighbors was 
first observed experi mentally and later understood 
microscopicallypl 10 ! 12 ! The main difference between the 
LDA-based models lies in the nearest-neighbor couplings 
J[ and J" that are sizable in the Br compound, yet rather 
weak in its CI counterpart. 

An important deficiency of the LDA-based approach is 
the neglect of FM contributions that are typical for short 
superexchange pathways. To evaluate the full exchange 
integrals (Ji), we turn to DFT+U calculations and list 



representative results in the last two columns of Table VI 
where we used certain values of the on-site Coulomb re- 
pulsion parameter U^, depending on the double-counting 
correction applied in DFT+U : Ud — 5 eV and 12 eV for 
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the around-mean-field (AMF) and fully-localized-limit 
(FLL) options, respectively. The use of a smaller Ud 
in AMF and a larger Ud in FLL is a rather general, al- 
though empirical, practice to ensure similar results ob- 
tained in the two approaches! 14 * 39 * ^ The specific choice 
of Ud is a more subtle issue, because the ratios and, par- 
ticularly, the energy scale of the couplings are strongly 
dependent on Ud- Following Ref.[HJ we adjust the Ud val- 
ues to experimental quantities, the Curie- Weiss tempera- 
ture 9 ~ 5 K that is a linear combination of all exchange 
couplings, and the saturation field HqH s ~ 70 — 85 T 
that measures AFM couplings in the system (see Ref. 6 
and Sec. IV C). However, the precise adjustment to both 



9 and H s with the single value of Ud is not possible, 
and uncertainties in the computed exchange couplings 
remain. The problem is probably unavoidable, owing to 
the complex crystal structure, large number of nonequiv- 
alent interactions, and their low energy scale. 

Combining the LDA-based evaluation of J^ FM and the 
calculation of Ji within DFT+J7, we arrive at a qualita- 
tive microscopic model of (CuBr)LaNb2 07. The main 
features of this model are: i) leading fourth-neighbor 
coupling J4 running through the particular, least curved 
Cu-Br-Br-Cu pathway; ii) large FM coupling between 
nearest neighbors along b: J1/J4 ~ —0.5; iii) sizable in- 
terlayer coupling J1/J4 — 0.3. Further details of the 
model remain ambiguous, and probably even obscure. 
The leading couplings J\ and J4 establish a stripe AFM 
order with parallel spins along b and antiparallel spins 
along a (see the right panel of Fig. 14). While weaker 



AFM couplings Ji, J", J 2 , and J4 support this ordering 
pattern, the possibly FM J' 2 and the small AFM J5 frus- 
trate the spin lattice. The relevance of these couplings is 
rather difficult to estimate, since their values depend on 
the computational procedure. For example, the sizable 
FM J' 2 appears in the FLL calculations only, while the 
AMF calculation results in a negligibly small J 2 . Unfor- 
tunately, the precise evaluation of such small couplings 
lies beyond the accuracy of present-day methods. There- 
fore, we further consider the problem of frustration using 
phenomenological arguments in Sees. |IVC] and [V] 



B. Magnetic structure 

Oba et alW have claimed the stripe AFM ordering 
in (CuBr)LaNb2 07, based on the observation of two 
magnetic reflections in a neutron diffraction experiment. 
They also reported the ordered magnetic moment of 
fi ~ 0.60 /is in surprisingly good agreement with the 
naive square-lattice model, which is, however, incon- 
sistent with the actual orthorhombic symmetry of a- 
(CuBr)LaNb2 07. We have revisited experimental infor- 
mation on the magnetic structure, because the powder 
data collected with the high-intensity D20 diffractometer 
allow to observe a larger number of magnetic reflections. 
Additionally, the refined value of [i depends on the scale 
factor for the nuclear scattering and, therefore, should be 
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FIG. 15. (Color online) Refinement of the magnetic structure 
for the subtracted neutron data (U.s k — I40 k): experimen- 
tal (circles), calculated (dark line), and difference (light line) 
patterns. Ticks show positions of the magnetic reflections. 
The excluded regions around 29 — 36 deg and 44 deg con- 
tain wiggles due to the incomplete subtraction of the 200/003 
and 202 nuclear reflections. Both subtracted and calculated 
patterns are offset from zero to ensure positive intensities. 



reconsidered for the revised structural model. 

Upon cooling below Tn — 32 K, at least two weak 
magnetic reflections were observed (Fig. |3|. To refine 
the magnetic structure, we subtracted the 40 K pat- 
tern measured right above T/v from the low-temperature 
pattern collected at 1.8 K. Since lattice parameters of 
(CuBr)LaNb2 07 are nearly unchanged below 40 K (see 
Fig. [7]) , the subtraction effectively eliminates the nuclear 
scattering and reveals a number of additional, very weak 
magnetic reflections at higher angles (Fig. 15). The sub- 



tracted pattern is refined as a purely magnetic phase. 

The magnetic reflections can be indexed by a k = 
(0,0, |) propagation vector that is consistent with k = 
(0, |, |) reported by Oba et aZPfor the four times smaller 
a su b x a su b x c su b tetragonal subcell. The Ah position of 
Cu and the k = (0, 0, |) allow for a number of irreducible 
representations, but only one of them leads to a complete 
refinement of the difference pattern. This representation 
corresponds to the stripe AFM structure, with the FM 
order of spins along either a or b. The spins follow the 
same direction, i.e., the a direction for the FM order 
along a and the b direction for the FM order along b. 
This magnetic structure is in agreement with the results 
reported by Oba et a/P The refined magnetic moment at 
1.8 K equals to 0.72(1) /j,b, which is slightly larger than 
the refined value of 0.60 /is reported in Ref. H3 



The stripe pattern (Fig. 15) features antiparallel spins 



along one of the intralayer direction (a or b) and along the 
intcrlayer direction c. The tiny difference between the a 
and b lattice parameters, as well as the limited resolution 
of the experiment do not allow us to determine the di- 
rection of stripes experimentally (Fig. [Id]). However, the 
microscopic magnetic model shows clearly that the cou- 
pling along the [CuBr] zigzag ribbons (b direction) is FM, 
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FIG. 16. (Color online) Possible magnetic structures with 
stripes of parallel spins directed along b (left panel) and a 
(right panel). Only the left structure is consistent with the 
NMR data, because Br atoms are strongly bonded to two 
Cu atoms with same direction of spin, thus featuring a large 
hyperfine field observed experimentally. 



while the coupling along a is AFM (TablejVT]). Therefore, 
we conclude that stripes of parallel spins run along the b 
direction. 

This conclusion is further confirmed by NMR results 
showing a large hyperfine field of 16.4 T at the Br siteP^ 
According to the structural models shown in Fig. [9j the 
Br atoms are surrounded by four Cu atoms, with shorter 
Cu-Br distances of about 2.5 A to the Cu atoms of the 
same [CuBr] ribbon, and longer distances of more than 
2.9 A to the Cu atoms of the neighboring ribbon. If 
stripes of parallel spins run along a, each Br atom will 
be strongly coupled to the spin-up and spin-down Cu 
atoms of the same ribbon, so that the hyperfine field on 
the Br site is strongly reduced or even canceled (Fig. 16 
right). By contrast, stripes of parallel spins running along 
the b direction (left panel of Fig. 16) ensure the strong 



hyperfine coupling of Br to the two closest Cu atoms 
with the same spin, and the ensuing large hyperfine field 
observed experimentally. 



C. Model simulations 

We will now consider experimental information on the 
magnetic behavior of (CuBr)LaNb20y and discuss the di- 
mensionality of the underlying spin lattice as well as its 
possible frustration. The DFT results (Table [Vl]) show 
clearly that the spin lattice is three-dimensional (3D), 
because sizable couplings are present in the ab plane and 
along the c direction (JjJ. The 3D nature of the spin 
lattice severely restricts the set of simulation techniques 
that could be applied to this system. Particularly, a fi- 
nite lattice of sufficient size can not be treated with exact 
diagonalization to perform a realistic simulation of a frus- 
trated 3D magnet. The method of choice is QMC that 
allows to handle large finite lattices, yet fails for frus- 
trated systems at low temperatures because of the sign 
problem. To perform QMC simulations, we neglect the 
potentially frustrating couplings J' 2 , J5, and J5, and con- 
sider the remaining interactions J\, J[, J", J 2 , J4, J4, and 
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FIG. 17. (Color online) Magnetic susceptibility of 

(CuBr)LaNb207 measured in the applied field of 0.5 T (cir- 
cles) and the fit with J 4 = 48 K, J1/J4 = -0.6, J2/J4 = 0.25, 
and J[ = J" = J2 = J4 = J± (solid line), see text for de- 
tails. The dashed line shows the intrinsic susceptibility of 
(CuBr)LaNb207 without the Curie- Weiss impurity contribu- 
tion Ci m p/ (T + dimp). The inset depicts an experimental mag- 
netization isotherm measured at 1.3 K (Ref.HH circles) and the 
fit with the same model parameters (solid line). 



J±. The reduction to the non-frustrated spin lattice is 
justified phcnomcnologically by the high Nccl tempera- 
ture (T/v/J 4 ~ 0.67, see below) and the sizable ordered 
moment of 0.72 (is, which is approaching the classical 
value of 1 hb for spin-^ (see the discussion in the end 
of this section). A posteriori, a remarkable agreement 
between the experimental data and simulations for the 
non-frustrated spin lattice also disfavors frustration in 
(CuBr)LaNb 2 7 . 

The spin lattice of (CuBr)LaNb207 features seven in- 
equivalent couplings that render fits to the experimen- 
tal data rather ambiguous. To reduce this ambiguity, 
we introduce the constraint J[ = J" = J 2 = J' 4 = J± 
since all these couplings are of the same order (see 



Table VI), and thus arrive at three variable parame- 
ters, J4, J\, and J2, as well as the g- value. Addi- 
tionally, we use: i) the temperature-independent term 
Xo to account for core diamagnetism and van Vleck 
paramagnetism contributing to the magnetic suscepti- 
bility; and ii) the Curie- Weiss term Ci mp /(T + 
to describe the low-temperature susceptibility upturrPH 
that presumably originates from weakly coupled impu- 
rity spinsP The temperature dependence of the suscep- 
tibility and the magnetization isotherm are both fitted 
with J 4 = 48(1) K, Ji/J 4 = -0.6(1), J2/J4 = 0.25(5), 
g = 2.29(2), xo = -7.4(3) x 10~ 4 emu/moL C imp = 
0.032(1) emu K/mol, and 6> imp = 6(1) K (Fig. fl7|) P Note 



that this set of exchange couplings perfectly reproduces 
the magnetic ordering temperature Tjv — 32 K, and the 
g-value is in the typical range for Cu +2 systems.^ Be- 
cause the g-value exceeds 2.0 and the saturated magne- 
tization is M s = gSfiB — 1-15 /is/f.u. for S = \, our 
estimate for the saturation field jJL$H s ~ 85 T is some- 
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what larger than /JoH s ~ 70 T reported in Ref. |H1 

The fitted value of Ci mp corresponds to about 8 wt.% 
of an impurity that shows non-negligible AFM interac- 
tions with the energy scale of #i mp = 6 K. Since our XRD 
and neutron data do not reveal any crystalline impuri- 
ties, the large Curie- Weiss term should be attributed to 
defects that are possibly related to random out-of-plane 
displacements of Br atoms. Presently, we have not stud- 
ied the effect of such displacements microscopically, be- 
cause even the model based on the equilibrium crystal 
structure is rather complex. The role of the out-of-plane 
displacements should be addressed in future studies. 

Despite the large number of variable parameters in the 
fitting, their estimates arise from different features of the 
experimental data, and appear to be robust under the 
constraint J[ = J" = J% = = J± introduced for 
individual exchange couplings. For example, C; mp and 
$imp are determined by the low-temperature part of the 
susceptibility that has a nearly temperature- independent 
intrinsic contribution of (CuBr)LaNb207 (see Fig. 17 1. 



The values of xo an d g rest upon the high-temperature 
data, where the susceptibility follows the Curie- Weiss 
law. Finally, the overall energy scale given by a linear 
combination of AFM couplings is fixed by the saturation 
field (slope of the magnetization curve), and the shape of 
the susceptibility above T/v determines J1/J4 and J2/ Ja- 

Once the constraint on the exchange couplings 
is released, the ambiguity appears. Similar to 
(CuCl)LaNb 2 07, the couplings can be largely redis- 
tributed between the weaker AFM exchanges J[, J", J 2 , 
and J4. Unfortunately, at the present stage this ambi- 
guity is unavoidable, given the large number of inequiv- 
alent exchange couplings and the limited experimental 
data collected on powder samples. A further refinement 
of individual parameters would require inelastic neutron 
scattering measurements on single crystals that provide 
access to the whole spectrum of magnetic excitations. 

We also evaluated the ordered magnetic moment and 
compared it to the experimental estimate (Sec. IV B I. 
Following Ref. 45] we calculated the staggered magne- 
tization (m s ) for finite lattices of different size up to 
24 x 24 x 12 and performed finite-size scalingP3 The 
m s values were obtained from the static structure factor 
taken at the propagation vector of the magnetic struc- 
ture, and from the spin-spin correlation taken at the 
largest separation on the finite lattice. The two ap- 
proaches yield the staggered magnetization of 0.879 \ib 
and 0.880 \ib and confirm the 3D nature of the spin sys- 
tem (compare to m s = 0.60 \xb and 0.83 \ib for the 
spin-i square lattice^ and cubic lattice,^ respectively) . 
To compare this result with the experimental fi deter- 
mined by neutron diffraction, one has to scale m s with 
the g- value (/i = gSm s ), and account for the Cu-Br 
hybridization that spreads the spin polarization to the 
ligand site, thus reducing the observed magnetic mo- 
ment on Cu. Using g — 2.29 and the DFT+U results 
showing the 20 — 35 % reduction in the moment due 
to the hybridization,^! W e estimate the ordered moment 



/i = 0.65 — 0.80 /is in remarkable agreement with the 
experimental result of 0.72(1) hb- 



V. DISCUSSION 

Our study provides detailed information on the crystal 
structure and magnetic behavior of (CuBr)LaNb207. We 
will now compare (CuBr)LaNb207 to (CuCl)LaNb207, 
and discuss their structural differences as well as the 
origin of different magnetic g round states. In agree- 
ment with earlier expectations we find that the CI 
and Br compounds have similar crystal structures fea- 
turing a tilting distortion within the [LaNb2 07] per- 
ovskite slabs and ordered displacements of Cu and halo- 
gen atoms in the ab plane. However, a salient feature 
of (CuBr)LaNb207 is the split Br position in the av- 
eraged crystal structure and an additional short-range 
order of Br resulting in the doubling of the c lattice pa- 
rameter. The Br atoms are also prone to out-of-plane 
displacements that become sizable at elevated tempera- 
tures (Fig. (8). 

The out-of-plane displacements represent deviations 
from the equilibrium crystal structure and can be largely 
understood as a dynamic effect, because the magnitude 
of the displacement increases with temperature (Fig. Is]). 
These displacements involve a small change in the to- 
tal energy of the system (Fig. 11 ), hence they are likely 
driven by soft phonon modes, although a comprehensive 
study of the phonon spectrum would be necessary to ver- 
ify this conjecture. The splitting of the Br position is an 
intrinsic feature of (CuBr)LaNb 2 07 driven by the small 
energy difference between two possible arrangements of 
the [CuBr] ribbons in the adjacent layers (Fig. [9]). While 
the coordination preference of Cu atoms ensures similar 
arrangement of all ribbons in the ab plane, the disorder 
between the adjacent planes has little influence on the 
structure and only slightly changes the energy. There- 
fore, diffraction data show the short-range order of Br 
evidenced by the diffuse scattering that reflects the pref- 
erential ABAB-type arrangement of the [CuBr] ribbons 
(Sec. [mil . 



The effect of split halogen position in the averaged 
crystal structure and the ensuing short-range struc- 
tural order may be relevant for all compounds of the 
(CuCl)LaNb 2 07 family. For example, single crystals of 
(CuCl)LaNb 2 07 also revealed two positions of CI atoms, 
although one of them was preferentially occupied (about 
90 % and 10 % of CI atoms, respectively) The increase 
in the c lattice parameter should further reduce the en- 
ergy difference between the AAAA and ABAB configura- 
tions, thus randomizing the arrangement of the [CuX] (X 
= CI, Br) zigzag ribbons in adjacent planes. This effec- 
tively eliminates the respective superstructure reflections 
and complicates experimental structural studies. For ex- 
ample, recent neutron diffraction experiments^ did not 
show any signatures of the superstructure formation in 
(CuBr)Sr2Nb30io, although the magnetism of this com- 
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pound, including the peculiar ^-magnetization plateaupSl 
can be hardly understood in terms of the square-lattice 
magnetic model derived from the parent, tetragonal crys- 
tal structure. A plausible explanation is the local order- 
ing of Cu and Br atoms that will not produce sharp su- 
perstructure reflections, which would be observable in an 
x-ray or neutron diffraction experiment. 

The effects of short-range structural order may also 
be relevant for solid solutions obtained by the CI/Br 
substitution. For example, CS2CUCI4 and CS2CUB1-4 
are orthorhombic compounds with fully ordered crys- 
tal structures, but the Cs2CuCl4_ :r Br a ; solid solutions 
have the tetragonal symmetry and presumably contain 
split positions of CI and Br at 1 < x < 2 (Ref. 2J. 
This tendency to the structural disorder should not be 
overlooked in the ongoing studies of quantum magnets, 
where the CI/Br substitution is used as a handy tool 
for changing the parameter regime and introducing bond 
randomness. The Br atoms not only occupy a crys- 
tallographic position different from CI, thus introduc- 
ing random exchange couplings, but also deteriorate the 
long-range structural order in the compound. Particu- 
larly, the (CuCli_ a; Br :c )LaNb207 solid solutions p2 show- 
ing the transition from the gapped magnetic phase at low 
x toward the long-range AFM order at high x, can not be 
considered as a simple experimental example of a quan- 
tum phase transition driven by a change in microscopic 
parameters. Structural changes are an equally important 
ingredient of the physics and require a further elaborate 
investigation. 

Despite the large spatial separation between the mag- 
netic [CuBr] units, (CuBr)LaNb2 07 does not show any 
experimental signatures of low-dimensional magnetic be- 
havior. Further, numerous couplings in the ab plane do 
not induce any notable frustration. Spin systems with 
reduced dimensionality, as well as frustrated spin sys- 
tems, feature strong quantum fluctuations that impede 
long-range magnetic order and trigger short-range order 
above Tjv- In antiferromagnetic systems, a typical signa- 
ture of the short-range order is the broad and symmetric 
susceptibility maximum, which is, for example, present 
in (CuCl)LaNb 2 7 and (CuCl)LaTa 2 7 (Ref. [55) yet ab- 
sent in (CuBr)LaNb2C>7 (see Fig. 17). Upon cooling, the 
susceptibility of (CuBr)LaNb207 drops down right below 
the magnetic transition at Tjy ~ 32 K, thus rendering the 
susceptibility maximum narrow and asymmetric. 

The short-range magnetic order evidenced by the sus- 
ceptibility maxima should be also observable in neutron- 
diffraction experiments above Tjy as a broad feature 
(diffuse scattering) preceding the formation of mag- 
netic reflections. The lack of diffuse magnetic scatter- 
ing at 40 K (see Fig. [3]) is another argument disfa- 
voring the low-dimensional and/or frustrated nature of 
(CuBr)LaNb2 07. Finally, the relatively high ordered 
magnetic moment fi ~ 0.72 [Lb is also a strong evidence of 
suppressed quantum fluctuations. In Sec. |IV C[ we have 
shown that the staggered magnetization m s ~ 0.88 /j,b 
(i.e., the ordered magnetic moment corrected for the hy- 



bridization with ligand orbitals and for the spin-orbit 
coupling) is similar to m s ~ 0.83 hb for the cubic spin 
lattice, the archetype 3D spin system with weak quantum 
fluctuations. 

In (CuBr)LaNb 2 07, weak quantum fluctuations 
should be ascribed to the large number of non-frustrated 
exchange couplings per magnetic site. Each Cu atom is 
involved in as many as ten couplings, compared to, e.g., 
six and four couplings on the cubic lattice and square 
lattice, respectively. Although quantum fluctuations re- 
cede because of the enhanced connectivity, they are not 
removed completely, as the reduction in the ordered mag- 
netic moment is partially related to the AFM exchange 
and the quantum behavior of Cu +2 (spin-i). For ex- 
ample, the staggered magnetization is reduced for 12 % 
compared to the classical value of 1 fig. Using classi- 
cal Monte Carlo simulations for the same spin lattice 
and same parameter regime, we arrive at Tn / J4 ~ 0.8 
and Tjv — 38 K compared to the experimental values 
of T N /J 4 ~ 0.67 and T N ~ 32 K. Therefore, weak 
quantum effects are still present in the spin system of 
(CuBr)LaNb207. This compound may be interesting as 
a system lacking both low dimensionality and strong frus- 
tration, so that weak quantum effects in a spin-i mag- 
net can be observed. In this respect, (CuBr)LaNb2C>7 
is rather unique, because Cu +2 compounds are prone to 
the formation of low-dimensional spin systems. In rare 
cases of geometrically 3D spin lattices, as in the green 
dioptase CugSigOis • 6H2O, the low coordination number 
(weak connectivity) still triggers strong quantum fluc- 
tuations that are clearly observed experimentally.^ The 
spin lattice of (CuBr)LaNb2 07 not only reveals the 3D 
geometry (the sizable coupling J± along the c direction) , 
but also features weak quantum fluctuations. 

(CuBr)LaNb 2 7 and (CuCl)LaNb 2 7 feature similar 
spin lattices, but strongly differ in the relevant param- 
eter regimes. The main ingredient of the spin lattice is 
the strong fourth-neighbor AFM coupling J4. This cou- 
pling runs via the long-range Cu-X-X-Cu superexchange 
pathway (X = CI, Br) and strongly depends on the devi- 
ation of this pathway from the straight line. Therefore, 
the second fourth-neighbor coupling is much weaker 
than J4, and spin dimers rather than bond- alternating 
spin chains are formed (see the right panel of Fig. 14). 



Although fits to the experimental data retain a certain 
ambiguity with respect to smaller individual exchange 
couplings, they yield a robust estimate of the effective 
interdimer coupling given by a linear combination of all 
interdimer couplings taken with appropriate coordina- 
tion numbers: J cS = 2| J x \ + J[ + J'{ + 2J 2 + J' A + 2J±. 
In (CuCl)LaNb207, J e g/Ji — 0.5 so that the system is 
close to the limit of isolated spin dimersP^ Therefore, 
(CuCl)LaNb2C>7 reveals a sizable spin gap^ and under- 
goes Bose- Einstein condensation of magnons in high mag- 
netic fields.^ On the contrary, (CuBr)LaNb2C>7 features 
Jeft/Ji — 3 and belongs to the opposite limit of strong 
interdimer couplings. The interdimer couplings close the 
spin gap and establish the long-range AFM order with a 
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relatively high T N /J 4 ~ 0.67. 

The strongly enhanced interdimer couplings are the 
crucial microscopic difference between (CuBr)LaNb2 07 
and its Cl-containing counterpart. Although DFT results 
are less accurate than fits to the experimental data, com- 
putational estimates show a similar trend of increased 
interdimer couplings in the Br compound and, particu- 
larly, reveal sizable nearest-neighbor interactions J[ and 
J" that are missing in (CuCl)LaNb2C>7 (compare Ta- 
ble VI to Table I in Ref. 14). The strong interdimer ex- 



change can be ascribed to the larger spatial extent of Br 
4p orbitals compared to CI 3p orbitals. The Br atoms en- 
able stronger interactions via short superexchange path- 
ways J{, J", and J2 that lack direct connections between 
the Cu0 2 X 2 magnetic plaquettes (see the left panel of 
Fig. 14). Moreover, the larger size of Br reduces the Cu- 
X-Cu angle within the zigzag ribbon (compare 101.8° to 
109.0° for Br and CI, respectively), thus increasing the 
FM coupling Jj. 

In summary, we studied the crystal structure, elec- 
tronic structure, and magnetism of (CuBr)LaNb207. 
This compound is very similar to its CI counterpart in 
general, yet different in several important aspects. First, 



the Br atoms are detrimental for the structural ordering 
and develop a short-range order, only. Second, the Br 
atoms are prone to out-of-plane displacements, especially 
at elevated temperatures. Third, the introduction of Br 
does not change the nature of the dimer-based spin lat- 
tice, although the interdimer couplings are increased dra- 
matically, thus triggering the transition from a gapped 
ground state in (CuCl)LaNb20y toward the long-range 
stripe AFM order in (CuBr)LaNb 2 7 . 
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TABLE S2. Details of structure refinements for (CuBr)LaNb 2 C>7. 
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FIG. SI. Rietveld refinement of the low-temperature D2B neutron data. The excluded region 
around 29 = 40 deg is due to the cryostat window. 
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FIG. S2. Rietveld refinement of the low-temperature D20 neutron data. The excluded regions 
around 29 — 18 deg and 40 deg are due to the magnetic scattering and diffuse scattering, 
respectively (see Fig. 6 of the manuscript). 
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FIG. S3. Rietveld refinement of the low-temperature synchrotron XRD data. The broad 
feature around 29 — 1.5 deg is due to the cryostat window. 
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FIG. S4. Rietveld refinement of the room-temperature E9 neutron data. 
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FIG. S5. Rietveld refinement of the room-temperature synchrotron XRD data. 
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FIG. S6. Rietveld refinement of the synchrotron XRD data collected at 720 K. 
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FIG. S7. Thermogravimetric data for (CuBr)LaNt>207. Note the onset of the weight loss 
around 750 K. 



